Proof Of The Invalidity Of The Bohzmann Property In The FMO Many-Body 

Neutrino Model 



James QuaclQ 

School of Physics, The University of Melbourne, Victoria 3010, Australia 

There has been recent debate over the use of the Boltzmann property in the kinetic equations 
describing dense neutrino systems such as early Universe and Supernova core. A technique de- 
veloped by Bell, Rawlinson, and Sawyer utilises the flavour evolution timescales of the neutrino 
systems to test the validity of this assumption. The Friedland-McKellar-Okuniewicz (FMO) many- 
body neutrino model was developed to conduct this test. It was concluded by its authors, using 
the Bell-Rawlinson-Sawyer timescale test, that the model lent support to the Boltzmann property 
- assumption. We developed kinetic equations for the FMO model. By direct analysis of the kinetic 

I equations we find, in stark contrast to Friedland et al., that in fact the Boltzmann property assump- 

. tion does breakdown in the FMO model. We have shown that the Bell-Rawlinson-Sawyer timescale 

' technique can only be used to invalidate the Boltzmann property but not validate it. 

m 

Single-body kinetic equations which encompass neutrino oscillation as well as neutrino-neutrino scattering have 
Qh' been found by Pantaleone [9l|,p^, Sigl and Raffelt [ll|, McKellar and Thomson 7] to describe neutrino dynamics in 
, early Universe and Supernovae. Underpinning these kinetic equations is the typical assumption, originally put forth 
by Boltzmann for classical dilute gases, that the neutrino density operator can be expressed as a tensor product of 
the individual one-body density operators of the system i.e. 
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I. INTRODUCTION 

This paper is based on the author's 2006 University of Melbourne BSc honour's report supervised by Bruce H.J. 
McKeUar. 
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, The validity of this assumption is not quantum mechanically obvious as interactions between neutrinos may produce 

■ entangled states which are not products of individual wavefunctions. 

' The effects of entanglement have up until now been looked from two approaches. The first by Friedland and 
. Lunardini [3(], constructed the flavour evolution of a many-body neutrino system as an effect of the interference of 
^ ' many elementary neutrino-neutrino scattering events. With this construction, they argued that in the limit of an 
. ^ infinite number of neutrinos the system can be factorised into one-body states. 

The alternate approach, introduced by Bell et al. studies the rate at which the many-body neutrino system 
H [ reaches statistical equilibrium. Essentially the argument here is that the one-body^ description predicts equilibration 

■ to be reached in time scales of an incoherent^ nature. Bell found equilibration to be reached in time scales of a coherent 
nature, thus citing evidence of the breakdown of the one-body description. Subsequently, Friedland and Lunardini 
Q found for systems with two types of neutrinos in equal number, equilibration time scales to be of an incoherent 
nature, providing support for the one-body description. Later, Friedland, McKellar, Okuniewicz generalised this 
model, and concluded equilibration time scales to be the incoherent type. This latter model is what we call the 
Friedland-McKellar-Okuniewicz many body neutrino model or FMO model for short. Because equilibration times 
scales were of an incoherent nature Friedland et. al. found, at least in the FMO model, no evidence for the failure of 
the one-body description. 

In this paper we enquire about the vailidity of the one-body description from a different angle. At the heart of 
the one-body description is the assumption that correlations between particles are insignificant. The question that 
naturally arises then is exactly what role does correlation affects have and are they indeed insiginificant. In an attempt 
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^ By one-body we mean factorisation into one-body states. This is an abbreviation used in the literature, which we will also use. 
^ Here, coherent and incoherent refers to the type of scattering 
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to answer this question we develop an exact kinetic equation of the Bogolubov-Born-Green-Kirkwood-Yvon (BBGKY) 
type. This wiU be used to measure the correctness of the FMO calculations. If they are correct, their results should 
exactly satisfy the exact kinetic equation. Then incorporating the one-body approximation, a Boltzmann-like kinetic 
equation will be developed. The degree to which this approximate equation is satisfied in the FMO model will then 
be investigated. The intention is that it will yield insights into the role of the correlation matrix and the validity of 
the one-body approximation in the FMO model. 



II. FRIEDLAND-MCKELLAR-OKUNIEWICZ MANY BODY NEUTRINO MODEL 



The FMO Model is a 'toy' model that was used by its authors to analyse the flavour evolution of a system of 
interacting neutrinos. So that the system could be exactly solved the authors made the simplifications that there 
are only two neutrino types and no spatial dependence. Although obviously this is not a representation of reality, it 
allowed the authors to exactly investigate the dynamics of flavour conversion. In this section we derive the hamiltonian 
and the one and two-body density matrices of the FMO many-body neutrino model. 



A. Hamiltonian Of The FMO Model 

The hamiltonian of the FMO model presented in this section is a reproduction of ^ albeit liberties with notation. 
In the FMO model, neutrino flavour is represented as spin. The interaction between neutrinos is described by the 
low-energy neutral current Hamiltonian 



Hnc — 




where = Val'^i^a 

a = is a flavour index 
G_F is the Fermi constant. 

Because this Hamiltonian is invariant under the rotation of the SU(2) flavour group the flavour space is isomorphic 
to spin interactions. 

Consider a two neutrino interaction. Since neutral current interactions conserve flavour, in the Hamiltonian of 
([2]) the flavour space wavefunction of a given outgoing neutrino is equal to that of one of the incoming neutrinos. 
Denoting the flavour wavefunctions as '^i and ^i, the only possible combinations are 



^*s,,^,nski^i (3) 

^*Su^inSk,^j (4) 
where we have suppressed the gamma matrices for clarity. 

Using the fact that X]i=i l^iP = S?=i I'^'iP — 1 the relation 2SiiSjk = SijSki + dlj ■ cfki, we can write the sum 
of (gD and © as 

1 + ^ (5) 

Thus the pair-wise interaction Hamiltonian can be written as. 



Hi. 
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J ^2 



(6) 



where g is a coupling constant. 

The interaction Hamiltonian for the whole system of N neutrinos is then, 
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where we have used the fact that 
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g[j^ + ^N{N-2) 



1=1 / i=l j=i+l 

B. The One And Two Body Density Matrices Of The FMO Model 

The original derivation of the one-body density matrix in the FMO model was conducted in [sj] . We reproduce the 
result here and extend it to the two-body case. 

1. Derivation Of The One Body Density Matrix 
We begin by considering the Hamiltonian of the system ([7]) ■ Its eigenvalues are 



H^j^M=g[j{J+l) + ^N{N-2)')^j,M (9) 

where V'j.a^ is an eigenstate of H. The subscripts remind us that the eigenstate is in the total angular momentum 
basis, J, and total angular momentum projection M . N is the total number of particles in the system. 
Let our system begin in the state 

\jumu > (^IjorriD >= [jujo^mumD > (10) 

where IJutt^u > is the total state of all spin-up particles and {Jd'tid > is the total state of all spin-down particles. 
Now the interaction Hamiltonian ([7]) will evolve this state in time. To do this conveniently, we wish to write our 
state in the total angular momentum basis, 

\jujD;mumD E < Jujo; J, M\mumD > iJujo; JM > (11) 
J 

Now we can easily use H to evolve this system in time, 



|*j,m(0 E < M\mumD;jujD > \juJd; JM > (12) 

J 

To get the one body density matrix, which we denote as pi, we need to single out this particular particle. To do 
this we decouple this one particle from the rest. 



Ijujo; JM >= \Uiku)ju,jD; JM > (13) 

where ji is the total angular momentum of the one particle, kjj is the total angular momentum of the rest of the 
spin-up particles. 
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Note that we have necessarily specified that the particle we are singling out is initially spin-up. An analogous 
derivation can be done for a particle that is initially spin-down. 

Now using known angular momentum coupling formula Eq. (jD5[l we write, 



\juJd; JM>^ hV^'"^'" [(2jc/ + l)(2j + 1)]^ 

.<j,j;m,m\JM>[f^ ^ ^^^^ 



X \jij;mim > 

So putting Eq. HI]) into Eq. ([H]) we have 

J.j,mi ,m 



X [(2jc/ + l)(2j + l)]^ 

X < iuio] JM\mijm]j >< jij; JM\mim > 



(15) 



Ji 3u 3n 
jo J 3 



\jij;mim > 



The density matrix of the total system is then, 



pit) = |*J,Af (0 >< 

J2 exp{-itg[J{J +1)- J' {J' + 1)]} 

m,m' ,mi ,m'^ 

X (2j„ + l)V(2j + l)(2j' + l) (16) 
X < jujo; JM\mumD >< jujn; J' M\mijmD > 

r 1/2 ju - 1/2 ju \ / 1/2 ju - 1/2 3U \ 
[3D J 3 j[ 3D J' 3' ] 
X \ jij;mim >< j[j';m[mi\ 

But we want the one-body density matrix. To do this, we take the trace over the rest of the system. 

pi(t)=Tr|j™>p(<) (17) 

Thus, 

P^' it) = E exp{-ztg[ J( J +1)-J'{J' + 1)]} 

m,mi ,m-|^ 

X (2j„ + l)(2j-f 1) 

X < jujD] JM\mumD >< jujn; J' M\mumD > (18) 
X < jij; JM\mim >< jij; j' M\m'im > 

r 1/2 Ju - 1/2 jN \ / 1/2 Ju - 1/2 JAT \ 
\3D J 3 ] \ 3D J' 3 ] 
X IjiTOi >< jim'il 

The superscript (f) in Eq. indicates that the particle denoted by the subscript, in this case particle 1, was 
initially in the spin-up state. For a particle that is initially down, one need only replace the U with D and vice versa, 
i.e., 
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p(^) it) = exp{-ztg[ J( J + 1) - J'iJ' + 1)]} 

m,mi ,m-L 

X (2j„ + l)(2j + l) 

X < joju; JM\mDmu >< joju] J'M\mDmu > (19) 
X < jij', JM\mim X jij; J' M\m'^m > 

^ r 1/2 jD - 1/2 jN 1 f 1/2 jD - 1/2 1 
\ iu J i J I Jc/ j J 

X IjiTOi >< jim'il 

^. Derivation Of The Two Body Density Matrix 

In this section we derive the two-body density matrix. If the two particles are initiaUy spin ahgned, the derivation 
follows that of i jllB 1) with the difference being that one must decouple two particles instead of one. The final result 
is that one need only make the following changes to Eq. (llSp . 

• replace ji with ji2 where ji2 is the summation of the angular momentum of the two particles i.e. ji2 = ji + j2 

• replace mi with toi2 where mi2 = -ji2, -ji2 + 1, ji2 - 1, ji2 



m,mi ,m']^ 

X (2i„ + l)(2i + l) 

X < jujo; JM\mumD >< jujo] J'M\mumD > (20) 
X < ji2i; JM\mi2m >< ji2j; J'M\m[m > 
1/2 ju - 1/2 JN 1 f 1/2 ju - 1/2 Jn 
in J J J I 3D J' j 
X |jl2™12 >< ji2m[\ 

P[2^ it) = E exp{-^^5[J( J + 1) - J'( J' + 1)]} 

m,mi ,m'j^ 

X (2i„ + l)(2i + l) 

X < joju; JM\mDmu >< joju] J' M\mDmu > (21) 
X < ji2i; JM\mi2m >< ji2j; J'M\m[m > 

1/2 JD - 1/2 JN \ ( 1/2 JD - 1/2 JN 
ju J j } \ ju J' j 

X |jl2TOl2 >< jl2m[\ 

The essence of the derivation for particles that are initially in opposite spins is to decouple one particle from the 
spin up system, and one from the spin-down system. To do this, we use the following three relations, 

Ijujo; JM >- E < ^ujd; JM\ji2j; JM > \j12j; JM > (22) 

ii2,i 



X 



\j12j, JM >= E ^ -^i^j; mumljuj; JM > \j12j; 17112111 > 



(23) 
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And with the help of Eq. jOS 



< {jiku)ju, {j2,kD)jD; JM\(jij3)ji2, {kuko); JM 

1 f Ji ku ju I (24) 

[ jl2 j J 

If we put the above three relations into Eq. ([12]), and recall that p[t) — |5'j,m(^) >< ^,/3f(i)|, we now have a 
density matrix in a form which is trivial to extract the two-body density matrix by tracing over the rest of the system. 
The final result is, 

p^2^\t)= exp{-^^5[J(J+l)- + 

J,^' ,j,ii2,i!.2 

m,mi2 ,?Tii2 

X {2ju + l){2jD + l)(2j + l)[(2ji2 + l){2j[2 + 1)]^ 

X < ju3D\ JM\mumD >< jujo; J'M\mumD > ^^g-j 
X < ji2j;rni2m\JM X jjjj; mi2m| J'Af > 

ji ku ju 
32 ko jo 

X \ji2m12 >< ji2»7ti2l 



Jl 


ku 


3u 


32 


ko 


3D 


312 


3 


J' 



C. Summary Of Parameter Definitions And Density Matrix Notation 



For convenience we summarise here the parameters used to describe the FMO model. 

• J is the total angular momentum of the system 

• M is the the projection of J 

• j is the total angular momentum of the rest of the system e.g. in the one-body density matrix it is the 
total angular momentum of the system less particle 1; for the two-body density matrix is is the total angular 
momentum of the system less particle 1 and particle 2. 

• TO is the angular momentum projection of j 

• ju is the total angular momentum of all the spin-up particles in the system 

• mu is the the projection of ju 

• j_D is the total angular momentum of all the spin-up particles in the system 

• TO£) is the the projection of jo 

• ku is the total angular momentum of the rest of the spin-up system e.g. in the one-body density matrix where 
particle 1 is initially spin- up, it is the total angular momentum of the system less particle 1; for the two-body 
density matrix of where particle 1 and 2 are initially spin-up, it is the total angular momentum of the spin-up 
system less particle 1 and particle 2. 

• ko is the total angular momentum of the rest of the spin-down system e.g. in the one-body density matrix 
where particle 1 is initially spin-down, it is the total angular momentum of the system less particle 1; for the 
two-body density matrix of where particle 1 and 2 are initially spin-down, it is the total angular momentum of 
the spin-down system less particle 1 and particle 2. 

• is the total number of particles in the system 

• (7 is a coupling constant 
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III. KINETIC EQUATION IN FMO MODEL 

Classically the time evolution of the phase space distribution is described by the Liouville equation, 

(26) 

where H is the Hamiltonian of the system, and curly brackets are Poisson brackets. 

The quantum mechanical analogue is the time evolution of the density matrix, which is described by the Von 
Neumann equation, 

'^-iH,p] (27) 
We can split the Hamiltonian into its single particle and interaction constituents, 

^ 1 

1=1 i=ij 

where K is the single particle Hamiltonian of the ith particle and Vij is the interaction Hamiltonian between the 
ith and jth particle. 

The rate of change of a single particle density matrix is found by tracing over all the other particles. 



= [Ki,Pi] + IJ2 [Vii , Pij] + J E [V^l , 



2 ^ --JL-j,r.jj . 2 
1 ' ' (29) 

9 E^^'yt^y'^iy] 



[if 1 ,Pi]+Y^ Tr, [V,, , pi,] + Tr^J- [^^J- ' P^^^^ 



2 

In the last line of Eq. (j29p we have employed the symmetry of our interaction Hamiltonian. 

Eq. (|29]) is in fact the quantum analogue of the first of the Bogolubov-Born-Green-Kirkwood-Yvon (BBGKY) 
equations @. 

In the FMO model there is no single particle Hamiltonian, so we will omit it. 

For the interaction Hamiltonian of our model it can readily be shown that the last term goes to zero. Here is the 
proof. 

Ti^jlVij , piij] = Try(< aia.iaj\Vtj piij\a[a'ia'j > - < aiaiaj\pujVij\a\a'^a'j >) 
= ^ < aiaittjlVij piijlu'ia^aj > - < aia.iaj\piijVij\a[aia'j > 



= ^ < aiaiaj\Vij\bibibj >< bibibj\piij\a[a'^a'j > 



(30) 



ai,aj 
bi .hi,bj 



- < aiaiaj\piij\bibibj >< bib^bjlVijla'^a'^a'j > 
where summation over primed equivalents is implied. 

With the appropriate choice of an orthonormal basis one can see that ai = bi = a'l, otherwise the term is zero. So, 
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Trij[Vij , puj] = ^ < aiaiaj\Vtj\aibibj X 016^6-, |aia-a^ > 

ZS- (31) 
- < aiaiaj\puj\aibibj >< aibibj\Vij\aia'ia'j > 

Using the fact that our interaction Hamiltonian is symmetric and that the appropriate basis consists of only two 
states for each particle, namely spin- up or spin-down, the summation over all possible ai, aj ,bi,bj , will yield zero as 
the negative sign in the second term will ensure that every term will cancel. 

Let's now consider the second term in Eq. i.e. ^ Trj[Vij, pij]. In this term there are two types of two-body 
density matrices: that where the spin of the two particles are initially aligned and that where the spin of the two 
particles are initially anti-aligned. Thus, Eq. can be simplified to. 



i-^ = NTr2[Vi2,Pi2]+MTr3[Vi3,pi3] (32) 

where we have labelled a particle that is initially anti-aligned to particle 1 as particle 2 and a particle that is initially 
aligned as particle 3. A'^ is the number of initially anti-aligned particles to particle 1, and M is the number of initially 
aligned-particles to particle 1. Note that N here is different from N in §11 C| which was defined as the total number 
of particles. 

In the next two sections we show that particles that are initially aligned to particle 1 do no contribute to its time 
evolution, and particles that are initially anti-aligned to particle 1 do. 

A. Initially Aligned Particles Do Not Contribute To Kinetic Equation 

The interaction Hamiltonian in matrix form is, 




(33) 



where U is the unit matrix. 

The spin basis of Eq. reading left to right (and of course top to bottom) is, | tT>5| ^i>^ iT>i| ll>^ where the 
first arrow represents the spin of particle one, and the second the spin of particle two. 
Let us now define a general two-body density matrix. 



Plr, 



' p21 p22 pl-i 
p31 pZ2 ^33 
,41 ^42 ^43 



n34 



Eq. ([331) has the same basis as Eq. Using Eq. and Eq, 

Trn\Vxri,Plr\ = 5 



we evaluate that, 



p32 _ p23 

„31 I „42 



^13 



,43 



P 



,34 
n32 



P 



p- - p"- + p- - p-" p- 

The two-body density matrix pi3 of Eq. (|32p in matrix element form is, 

< mi,TO3|pi3(t)|mi,m3 > 

= E (-)^-'''exp{-ztg[J(J+l)- J'(J'-Hl)]} 

J. j' ..j,7n 

X (2jt/ + l)(2j + l) 

X < jujo; JM\mumD >< jujo; J'M\mumD > 
X < jisj; JM\mi3m >< jisj; J' M\m\^m > 



Ji3 ]u - Jl3 Ju 
3D J j 



Ji3 Ju -313 Ju 
3D J' 3 



(34) 



(35) 



(36) 
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Now mi3 + m = M. So for element < ^Ipial^, —5 >, 

< jisj; mi3m\JM X jisj; m^^mlJ'M > 

Ij; lm\JM >< Ij; Om\j'M > (37) 

= 

since it is not possible for M to satisfy both mi3 + m = 1 + m and m']^3 + m = + m. 

Therefore matrix element =< i, ^Ipial^,— 5 >= 0. By same reasoning matrix elements, p^^ — p^^ = p^'^ ~ 
p3i _ p2i _ ^23 _ pi2 _ Q Yrom Eq. p5|) we thus see that the off-diagonal elements of A/Tr3[Vi3, P13] are zero. 

Let us now make the following definitions, 

a(J, J') ^ ^(-)'^-^'(2j^ + l)(2j + 1) (38) 

X < jcio; JM\mumD >< juiv, J' M\mumD > (39) 
X < jid; JM\mi3m >< ji3j; j'M|m'i3m > (40) 

^ f ii3 - ji3 ju } f ji3 - ji3 ju 1 /^-^^j 

b{J) = -itgJ{J+l) (42) 
Then pi3 can be written in the following way. 



pi3(0 = 51 "(^' J>MKJ) - K^')} 
./,,/' 

= ( a{J, J') exp{6( J) - b{J')} + a{J', J) exp{5( J') - 6(J)} j + ^ a( J, J' = J) 
\j<J' J J 

= '^'){cos(6( J) - b{J')) - i sin(6( J) - 6( J'))} + a(^', J){cos(&( J') - b{J)) - i sin(6( J') - b{J))}j 

+ Y.^{j,j' = j) 

(43) 

Now a( J, J') is symmetric about J, J' i.e. a( J, J') — a{J' , J) . Therefore, 

pi3(t) = ^ 2a( J, j')cos{b{J') - b{J)) + a( J, J' = J) (44) 



J,J' J 



Consider matrix elements < 5,— 5IP13I — 5, 5 > and < —5,^1^13!^,—^ >. Getting from the former element to 
the latter element is equivalent to swapping primed and unprimed variables. But from Eq. (|44p we see that pisit) is 
invariant under primed and unprimed variable swapping; hence < ^|/9i3| — 5,5 >=< 5IP13I5, ~^ >• This 
tells us the diagonal terms in Eq. p5|) are zero. 

Thus we have shown that the aligned term in the kinetic equation, AfTr3[Vi3, P13], is zero. 



B. Initially Anti- Aligned Particles Contribute To Kinetic Equation 

The argument for initially aligned particles that showed that p^^ = p^^ = p^'^ = = p^^ — p^^ = p^^ = = 
also hold for the case of initially anti-aligned particles. However careful consideration of Eq. pS)) reveals that it is not 
invariant under the swapping of primed and unprimed variables, therefore the two-body density matrices for initially 
anti-aligned particles is not symmetric. From Eq. (|35p we see that a non-symmetric density matrix will populate the 
diagonal terms, thus the anti-aligned term of the kinetic equation will contribute to the kinetic equation. 

Therefore the one-body kinetic equation reduces to, 

i^=NTr2[Vi2,pi2] (45) 

where N is the number of particles in opposite spin to particle 1, and particle 2 labels a particle that is initially in 
opposite spin to particle 1. 

From now on, we will just refer to the one-body kinetic equation as the kinetic equation. 
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C. Verification Of Consistency Of FMO Model With The Kinetic Equation 

As the kinetic equation is an exact equation, it stands to reason that the density matrices of the FMO model should 
satisfy the kinetic equation. In fact this would be important verification of the correctness of the one and two body 
density matrices of the FMO model. In this section we outline the proof that was used to show the consistency of the 
FMO density matrices with the kinetic equation. 

The proof consists of two parts. In the first part, it is shown that the left and right-hand-side are of the same form, 

.dpi (A 



iVTr„[yi„,pi„]= (47) 

In the second part through general angular momentum properties and physical constraints specific to the FMO model, 
it is shown that A — B. 

The complete proof of the consistency of the FMO density matrices with the kinetic equation can be found in 
Appendix [A] 



IV. ANALYSIS OF ONE-BODY FACTORISATION ON KINETIC EQUATION IN FMO MODEL 

Our investigation into approximating the density matrix as a tensor product of one-body density matrices begins 
by seeing the effects it has on the kinetic equation in the FMO model. The kinetic equation with factorised density 
matrix becomes, 

i^^NTr2[Vi2.Pi®P2] (48) 

Eq. (|48| is in fact analogous to the classical Boltzmann equation which is the transport equation for dilute gases. 
Implicit in the derivation of the Boltzmann equation is the Boltzmann property, 

N 

P=X{Pn (49) 
n=l 

the quantum analogue of which we are investigating. For brevity Eq. ()48p will be referred to as the Boltzmann 
equation in the rest of this paper. 

Let us start by revisiting the FMO one-body density matrix, which we present here in matrix element form, 

< rni\pi\m'i > 

= 9 exp{-%(J(J+l)- J'(J' + 1))} 

J, J',j,m 

X (2jy + l)(2j + l) ^^^^ 
X < iujo; JM\mumD >< jujo] J'M\mumo > 
X < Jii; JM\mim X jij; J'M\m[m > 

" 1/2 ju - 1/2 ju 1 f 1/2 ju - 1/2 ju 
in J J J I JD J' i 

Now the term < jij; JMjmim >< jij; J' M\m[m > forces the condition M = mi + m = m[ + m, which implies 
that the only density matrix element that is non-zero is when mi = m'l i.e. the FMO one-body density matrix is 
diagonal. Therefore the tensor product of the one-body density matrices i.e pi (g) pn is also diagonal. 

Now if we recall Eq. (I35p ). which says that for a general two-body density matrix, 

Trn[Vin, Pin] = 5 ( ^21 _ ^31 ^ ^42 _ ^43 ^23 _ ^32 j (51) 
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then, 

NTt2[Vi2,Pi®P2]^0 (52) 

as pi (E) P2 is a diagonal matrix. 

This was a result we did not expect. Of course in general is not zero. So what does this all mean? 

To answer this question we need to realise that the density matrix can be written as the sum of the tensor product 
of the one-body density matrix and the correlation matrix^ 



P ■ 



Yl^Pn + r (53) 



So if we rewrite the kinetic equation with the correlation matrix explicitly included, 

«^ = A^Tr„[Vi„, Pin] 

= NTvniVir^, Pi ® Pn] + ^Tr„[Fi„, r] (54) 

= iVTv„[yi„,r] 

we clearly see that the evolution of the system in the FMO model is entirely driven by the correlation matrix. 

Let's now take a moment and recall the works of Bell et. al [l| and Friedland et al.[3|[l]- The premise in these 
works was that the probability of flavour change due to coherent scattering is. 



= i?2|A|2 (55) 



and the probability of flavour change due to incoherent scattering is, 

i 

where A is the amplitude of a scattering event occurring and R is the number of background neutrinos. 
The time scales on which the probability evolves is, 

tcoh oc (57) 

tine OC (58) 

for coherent and incoherent scattering respectively. See for further details. 

Hence if the system exhibited flavour conversion in timescales which was not of an incoherent nature, then this was 
an indication of the breakdown of the one-body approximation, as the one-body approximation predicts conversion in 
incoherent timescales. Friedland found that in the FMO model, flavour evolution was of an incoherent nature, citing 
evidence of no breakdown in the one-body description. However what we have found here is that the flavour evolution 
of the FMO model is entirely driven by the correlation matrix, and in fact the one-body description would describe 
no flavour evolution at all; clearly there is a breakdown in the one-body description in the FMO model. 



A. Effects Of Adding Neutrino Vacuum Oscillations 



The standard FMO model only encompasses flavour evolution through neutrino-neutrino interaction. We wish to 
extend this model to include the neutrino vacuum oscillation. To achieve this, as originally suggested by Okuniewicz 
0, we introduce a constant magnetic field to the system. This works because we are representing neutrino flavour as 
spin, and the introduction of the magnetic field will cause these spin vectors to precess about the magnetic field. As 
the spin vector processes, the spin projection on to some axis which represents a flavour state will oscillate (as depicted 
in Fig[T]). This means the probability of being in some flavour state will oscillate representing the phenomenon of 
neutrino vacuum oscillations. The reason we are interested in adding neutrino vacuum oscillation to the model is that 
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FIG. 1: _B is the magnetic field vector about which spin vector P precesses. The spin projection on to the flavour axis determines 
the probability of finding a neutrino in that fiavour. The value of this probability will oscillate between a and b, the maximum 
and minimum value of the projection of P onto the flavour axis. 



it will populate the ofF-diagonal terms in the one-body density matrices - remember it was because of the diagonality 
of the one-body density matrices in the FMO model that had lead the factorised term to zero. 

The derivation of the density matrices with oscillations'^ follows in a similar fashion to the system without oscilla- 
tions, except that now the Hamiltonian has an extra single particle term, 

H = K, + (59) 

= B-J + gJ^ +g^N{N -2) (60) 

where B is the magnetic field vector, designed to be constant. The one and two body density matrices with oscillation 
can be found in Appendix [Cl 

Including oscillations, the kinetic equation is now 

i^^Kpi+NTr2[Vi2,Pi®P2] (61) 

With density matrices which are now in general not diagonal, the argument presented at the beginning of this section 
that shows that Tr2[Vi2, pi (8) P2] is zero no longer holds. Surprisingly however, a numerical analysis will maintain 
that this term is zero. This tells us that, at least for the numerical cases we have tried, that irrespective of whether 
we include vacuum oscillation in our model, the factorised term is still zero. The next section will show that in fact, 
in general, extending the FMO model to include vacuum oscillation will not alter the actuality that the dynamics of 
the system is completely driven by correlation affects. 



V. WHY THE ONE-BODY FACTORISATION BREAKS DOWN IN THE FMO MODEL 



To understand why the one-body description completely breaks down in the FMO model, we ask the converse 
question: when will the one-body description not completely breakdown? 
To ascertain this, let's rewrite the factorised term in the following way, 



In the literature neutrino oscillations implicitly mean neutrino vacuum oscillation. This is the terminology we will use also from now 
on. 
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(62) 



Let 



Then, 



Tr2[Vi2,pi (8)/32] = ^ < mim2\[Vi2, pi ® P2\m'lm2 > 

= ^ < miTO2|^12Pl ® P2- Pi® P2Vl2\m'(m2 > 

= ^ < mim2\Vi2\rn[rn'2 >< ml'\pi\m'( >< rn'2\p2\m2 > 

m2 .rn'^ ,in'^ 

— < mi\pi\m[ >< TO2|p2|'7i2 >< m[m'2\Vi2\m'(m2 > 

= ^ < mim2\Vi2\m[rn'2 >< rn'2\p2\'rn2 >< ml'\pi\rn'( > 

m2 .rn'^ ,mfj 

— < mi\pi\m[ X ■m\m'2\Vi2\'m!lm,2 >< m2\p2\m'2 > 



< mi\fi\m[ >= ^ < m[m'2\Vi2\m'(m2 X m2\p2\m'2 > (63) 

7712 -77l2 



(64) 



Tr2[Vi2, P1P2] = ^ < milfilm'^ X m[\pi\m'{ > - < mi\pi\m[ >< mi|/i|m'/ > 
=< mi\[fi,pi]\m[ > 

Now recall from mi Al that the interaction Hamiltonian between two particles is |(3 + iXi • cr^ ). 
So if we write the density matrix as a linear combination of the Pauli matrix basis, 

Pi = liPo + P ■ 'T,) (65) 

P2 ^ ^iQo + Q ■ (fi) (66) 

and realise that the only term in the interaction Hamiltonian that doesn't necessary vanish in commutation is ^cTi -o^, 
then fi can be expressed as, 



fi = Tr2{Vi2P2) 

= |rr2((Ti;(T2, 

= ^Tr2{(TiiC72iQo + criiQjniShn + «e/77ma'2n)) (-gy-j 



-Tr2{cril(72l{Qo + Qm(T2,n)) 



Hence, 



= |rr2(fTl • Q) 

9 ^ A 



Tr2[Vi2,Pi ® P2] = |[<? • g, Po + <y^i ■ P] 



= ^QlPmWlhf^lm] 

_ % (68) 
= i|(QxF).al 

Eq. (j68[) tells us that the evolution of the system will be driven entirely by the correlation matrix when the interacting 
spins are parallel, and conversely the factorised term will contribute to the dynamics of the system when the interacting 
spins are not parallel. Now we have found that for our model that the factorised term is zero for all time; this suggests 
that the interacting spins are always parallel. 
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Now the left-hand-side of kinetic equation can be written in the following way, 



5pi IdP (gg) 



Combining Eq. ([S5|) .Eq. the Boltzmann equation becomes, 

dP - - 

— = Ng[Q X P) (70) 

Eq. (jTO]) is of a familiar form; it is indeed the equation for precession. It tells us that that the time evolution 
of a self-interacting neutrino system can be exactly described by the factorisation into one-body states only if the 
interacting spins precess around one another. As the one-body description yields absolutely no dynamics, this clearly 
is not happening. Why this is so, is evident if we recall that our one-body density matrix in the FMO model is diagonal 
in the spin basis. What this means is that if you were to single out an individual particle, which you necessary do in 
the one-body description, you will only ever find it in the spin-up or spin-down state. Particles which are spin-up or 
spin-down are parallel and therefore their spin vectors can never precess around one another. 



A. Breakdown Of Oscillation Extended FMO Model 



The above analysis also explains why the factorised term in the Oscillation Extended FMO model is zero. Initially 
in the FMO model the system starts in an untangled state with every particle in either an up or down spin state. At 
i = when the magnetic field is switched on, spin-up particles begin to precess, say, clock-wise around the magnetic 
field vector B. Spin-down particles precess in an opposite direction, that is, anti-clockwise. We define (anti-)clockwise 
by looking along the direction of the vector that we are describing. As illustrated in Fig. [2] the particles spin vectors 
will remain parallel drawing two cones; hence the introduction of the magnetic field, although causing the particles 
to precess does not distort them from their parallel configuration. Therefore by Eq. ([70)1 the one-body description in 
the Oscillation Extended FMO Model also breaks down. 



VI. HOW THE KINETIC EQUATION EXPLAINS FLAVOUR EVOLUTION IN FMO MODEL 

We extend the analysis in Chapter |V] to show how the kinetic equation can represent the complex flavour evolution 
in the FMO Model found in the work of Fricdland et. al and OkuniewiczQ. 

The interaction Hamiltonian of the FMO model is a special case of a more general form, 

Vl2 = VijGuG2j (71) 

where j, j — 0..3, and summation over repeated indices are implied. 
We can also express the two-body density matrix in a similar fashion, 

P\2= Vi3^\i^2j (72) 

Now, 

= VijPkiicriiCrikcr2jCr2i — criiaikcr2jcr2i (73) 

+ criiCrifc(T2iCr2; -|- cr ikC ii(T2lcr2j) 
= Vi-jPki{<yiiCrik[a2j,(T2i] + [(Tii,crik]a2icr2j) 

Since aio and CT20 are just the identity operators, they disappear in the commutation. Thus, 

[Vi2,Pi2] = 2ivijPki{ejhnO-iiCrikcr2m + iikn<yin<y2io2j) (74) 
where i. A: > in the first term, and j,l > in the second term. 



15 



FlavDur Axis'- 



■* 
■1 



FIG. 2: Spin vector P began its life aligned with the flavour axis. Q began anti-aligned with the flavour axis. When B was 
switched on, P began to precess. Looking along the direction P, it began precessing clockwise. Looking along the direction of 
Q, Q also precesses, but in an anti-clockwise direction such that P and Q draw out two cones. 



Further, 



Tr2[yi2,pi2] = 2TOjjp«Tr2(ejiTO0-H(Tucr2m + eife„cri„0-2;cr2j) 

= ^iVijPkl<^iknCrin5lj (75) 
= 4:iViiPkieikn<7ln 



where i,k > 0. 

Before continuing let us make the following definitions, 



V = matrixiwy} (76) 

p = matrixjpij} (77) 

w = matrixjwfe;} (78) 

q = matrix{4pfei} (79) 



where i,j = 0..3 and kj — 1..3. 

The usefulness of these definitions will be apparent soon. 
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Now we know something about p. Firstly, 

pi = Tr2pi2 = 2pioaii (80) 
P2 = Tripi2 = 2pojCT2i (81) 

So expressing pi and p2 as a linear combination of the Pauli matrix basis as in Eq. (|65p and Eq. (j66p and applying 
the condition that the trace of density matrices equal 1, then we know that p must at least be of the form, 



P=4 



n Qx Qi Qs^ 

Pi Pll Pl2 Pl3 

P2 P2I P22 P23 

P31 P32 P33/ 



(82) 



It is prudent to be clear of the distinguishment between Eq. (|82|) and Eq. (|34p . The latter is a density matrix, whereas 
the former is a co-efficient matrix. Their exact relation is described in Eq. ([7^ . where Eq. (IM)) is the Ihs and Eq. ([5^ 
is the matrix co-efficient on the rhs. 

With this information about p in hand, we can express Eq. (|75p as, 

Tr2[yi2,Pl2] = ^iViOPkO^iknC^ln + 4:iVijPkjeiknCrin 

= i{v"o X P) ■ ai+ i{wq^)ikeikn(yin (83) 
= i{vo X P) • CTi + i{wq^)^ ■ ai 

where i, j, k,n — 1..3, and 

^'O = (^'I0,'f20,'f30) 

In Eq. ((83|l we have borrowed from exterior differential calculus the notion of a dual of a matrix i.e. for a general 
3x3 matrix M, the vector = djkMij.^ 

Therefore the kinetic equation can be written as, 

1 f)P 

^^ = (^roxP) + (,«g^)^ (84) 



In the FMO model. 



hence, 



V12 = ^3 + ^<^i<^2 (85) 



/3 0\ 
^FMO^ff 100 
^ 2 1 

Vo 1/ 

Again we make clear the difference between the two matrices Eq. and Eq. (IHS)) . The latter is the interaction 
Hamiltonian in matrix form, the former is a matrix co-efhcient. Their exact relation is described in Eq. (|7ip . where 
Eq. p3)) is the Ihs and Eq. ((86| is the matrix co-efficient on the rhs. 

From Eq. ([55]) we see that Vq = and w |ll. Therefore if we revisit Eq. ((84|) with the FMO Hamiltonian, the 
first term disappears but we are left with the second term, 

^ = Ng{q^r (87) 

The probability of finding a neutrino in a certain flavour state is determined by the projection of P on to the axis 
that represents that flavour state. Hence Eq. ((87)) is exactly what drives the the flavour evolution in the FMO Model. 



^ In mathematical text the dual is usually represented with *, e.g. M^. We have chosen to use ^ instead to avoid confusion with the 
complex conjugate which is also sometimes represented with a *. 
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Eq. ([87|) is the exact kinetic equation of the FMO model. As a check, we should be able to retrieve the factorised 
kinetic equation by applying the Boltzmann property, 

Pl2 = Pi® P2 

_lpn ^^^^ 

where i,j — 0..3 

Comparing Eq. (I88p with Eq. (17^ we see that pij = jPiQj. Therefore, 

= tlmnPmQl (89) 

= QxP 

where l,m,n = 1..3 

Substitution of Eq. ([SQ]) into Eq. ((87|) recovers exactly the factorised kinetic equation of Eq. ([70|) . 
Furthermore, from Eq. ([87|l we able to reveal that the FMO is driven entirely by the correlation term. To see this 
we first expand out {q'^)^, 



(90) 



I T T T T T T \ 
= [123 - 932 > 931 - 9l3> 1l2 ' 921 ) 

Now recall, 

Pi2 = Pi P2 + r (91) 
If we define the correlation matrix in the following way, 

r=^7yCrHCT2j (92) 

then Eq. (pijl can be written as, 

P^j = -PiQj + (93) 



4 — ■ 4 

With the constraints imposed on 7^ by Eq. (|82|) . we see that 



1 

^=4 



/I Qi Q2 Qs 

Pi PiQi + -fii P1Q2 + 112 P1Q2 + J13 

P2 P2Q1 + 721 P2Q2 + 722 -P2Q2 + 723 

V^'3 ^'3^1 + 731 P3Q2 + 732 ^'3^2 + 733; 



(94) 



Hence, 

{q^)^ = (P3Q2 + 732 - P2Q3 + 723, ^'ll93 + 713 - ^'3Ql + l31,P2Ql + 721 - P1Q2 + I12) , , 

Now we know that in the FMO model Q x P = 0, thus, 

% = A^.9(7^)^ (96) 
Therefore in Eq. (|96p we have shown how the flavour evolution in the FMO model is entirely driven by the correlations. 



VII. CONCLUSIONS 



The ability of physics to abstract complex systems to simple ones by stripping away that which is not essential has 
been extremely successful; it has allowed mankind to comprehend his Universe to levels of detail only dreamt about 
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in the past. However the strength of these simpHfied models ironically rests on the validity of what is not included 
in the model. Time and time again in physics' past, questioning the legitimacy of these assumptions have opened a 
Pandora's box of insights and new physics. 

To describe the kinetic equations of dense neutrino systems such as early Universe and Supernovae, the literature 
has assumed the Boltzmann property where many-body states can be portrayed as the product of the individual 
states of the system. Debate over the validity of this assumption in such systems has recently surfaced. A technique 
which has been used to test the legitimacy of this assumption is to compare the neutrino flavour evolution of the exact 
many-body systems to that which is predicted by the product of one-body states. Since the one-body description 
predicts flavour evolution to be of an incoherent timescale, if the exact system was found to have flavour evolution 
which was not of an incoherent nature, than this would signal the breakdown in the one-body approximation. The 
Fricdland-McKellar-Okuniewicz model was developed to conduct this test. It was found by the authors of the model 
that their system evolved in incoherent timescales, and hence they concluded support for the one-body approximation. 
Our investigations found that the dynamics of the FMO model was entirely driven by the correlation function, and 
in fact the one-body description would describe no evolution in the system at all; there clearly was a breakdown in 
the one-body description in the FMO model. 

It is important to recognise that the Boltzmann factorisation of the two-body density matrix was not the only 
assumption used in the kinetic equations of the early Universe and Supernovae. Another assumption central to the 
derivation of the neutrino kinetic equations was that the duration of a collision was short compared to the time 
between collisions. The same assumption was made in Boltzmann's derivation of his equation for the kinetics of 
dilute gases. In the FMO model there is not spatial dependence in the interactions, so that the interactions are never 
turned off. The duration of the interactions is thus the whole time elapsed. With this in mind, perhaps we should 
not be so surprised that the Boltzmann-type of kinetic equation breaks down in the FMO model. Notwithstanding 
the increased complexity, a better model to test the Boltzmann property with will need to include or at least mimic 
a spatial extent. 

Our conclusion that the one-body description breaks down in the FMO model is in stark contrast to the result 
found by Friedland et. al. This questions the effectiveness of using neutrino flavour timescales to test the validity 
of the one-body approximation. Even though the FMO model evolved in the same timescales as that predicted by 
the one-body approximation, this is not an indication of the validity of the approximation, as we have shown that 
this approximation in the FMO model would catastrophically breakdown. We conclude that at best, the technique of 
comparing neutrino flavour evolutionary timescales can only invalidate the Boltzmann property but never validate it. 
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APPENDIX A: PROOF OF THE CONSISTENCY OF THE DENSITY MATRICES WITH THE KINETIC 

EQUATION 

In this appendix we show that the one-body density matrix and two-body density matrix derived in i)II Bl is consistent 
with the kinetic equation Eq. (|^5)) . The proof is divided up into logical sections IbIIB 21 

APPENDIX B: FORM OF LEFT-HAND-SIDE 

Let's begin by differentiating Eq. to get the left-hand-side of the kinetic equation (note although we are 
considering the case when the particle's initial configuration is spin-up, it also holds for the case when the particle is 
spin down). 

'^=9 {-f''''\j{J + 1) - J' {J' + l))/exp{-itg{J{J + 1) - J'( J' + im23u + l)(2j + 1) 

J, J' j',m 

< iuio; JM\mumD >< jujo] J' M\mumD >< jij; JM\m\°'^m >< jij; J'M\ml°'^m > (Bl) 

r 1/2 ju - 1/2 ju \ / 1/2 ju - 1/2 ju \ 
\jD J J j \ JD J' J i 
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Here we have expressed the Ihs in matrix element form, m™™ and m™ represent the spin projection basis of the row 
and column of the matrix. One sees that the condition M = m™™ + m = m™' + m in Eq. (|B1[) implies that the Ihs 
is diagonal in this basis. 

Next we acknowledge the following two properties: 

• that for J ~ J', the terms in the summation go to zero since J(J + 1) — J'(J' + 1) = 

• from the triangular constraints in the 6j-symbol, J, J' G {|j — j + ^} 

If we combine these properties with the known Clebsch-Gordan co-efficient relation, 

< jij2mim2\jij2JM {~y^-3^-3^ < j2jim2mi\j2jiJM > (B2) 

we see that i-^-^ — = ' '^-'^si'S the superscript is the (m™™, m^^"') of Eq. (jBip . which identifies the matrix 

element. 

Let's now have a look at the rhs of the kinetic equation. 



1. Form Of Right-Hand-Side 



The interaction Hamiltonian in matrix form is, 



/l ON 
0-120 
2-10 

\0 1/ 



(B3) 



where U is the unit matrix. The spin basis of Eq. (|B3p reading left to right (and of course top to bottom) is, 
I TT>J T-i>il iT>i| i-l>7 where the first arrow represents the spin of particle one, and the second the spin of particle 
two. 

Let us now define a general two-body density matrix. 



Pl2 



/p" ^13 ^14. 

p21 p22 ^23 ^24* 

^31 ^32 ^33 ^34 

Vp41 p42 ^43 ^44 



(B4) 



where we have omitted the subscripts in the matrix elements for clarity. Eq. (|B4p has the same basis as Eq. (jB3 
Using Eq. (jBSp and Eq. (|B4p to evaluate the rhs of the kinetic equation we get. 



A^TV2[yi2,Pl2]=A^5(^21_^31 



p23 pl3 _ pl2 + ^34 

p42 _ ^43 p23 _ ^32 



P 



Now mi2 



= M. So for clement < ^, 5IP12I 



2 ' 2 



< iini\minm\JM >< j[.J;m[^m\J'M > =< Ij; lm\JM >< lj;Om\J'M > 

= 



(B5) 



(B6) 



since it is not possible for M to satisfy both mi2 + m = 1 + m and m'12 + m, = + m. 
By same reasoning matrix elements, p^^ = = p^'^ = p^'^ = p^^ — p^^ — p^^ — p**^ = 0. 

Using this result, the off-diagonal terms of the rhs matrix Eq. (|B5[) are zero. As for the diagonal elements, Eq. (|B5[) 
clearly tell us they are the negative of each other. 

To this point we have shown that the Ihs and rhs of the kinetic equation exhibit the same form, namely the 
off-diagonal matrix elements are zero, and the two diagonal elements are the negative of each other. The only thing 
left to do now is to show that one of the diagonal elements of the Ihs is equal to the equivalent one of the rhs. 



1. Showing The Equivalence Of Ihs and rhs Diagonal Elements 

We are going to show that Ihs^^'^^ = rhs^^'^\ where the superscripts have their usual meaning, in this case 
identifying the first element of the matrices. 
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3. Ihs^i'i^ 

Define, 

a{J,J') = JiJ + l)-J'{J' + l) (B7) 



/3(J,J') = (2jc/ + l)(2j + l) 

, , 1 ,1 1 , ,1 

< jujo; JM\mumD >< jujo] J M\mumD >< -j; JM\-m >< -j; J M\-m > 

r 1/2 ju - 1/2 ju 1 f 1/2 ju - 1/2 ju 

\ JD J j i \ JD J' j 

Then, 

Z/is(5'5) ^ a{J,j'){cos{tga{J,j'))-isin{tga{J,J')))(3{J,J') 

= -2ia{J,J')sin(tga{J,J')))(3{J,J') 

J<J', 
j 



(B8) 



(B9) 



The 6j-symbol in (3 restricts J,J' € {\j — 5I, J + |}; and because J < J',then J = |j — = J + ^- Note that if 

J = J' then a = 0. With these substitution for J, J' we get, 

Ihs^--^--^ = i2g{2ju + 1) £ (2j + 1)2 sm{tg{2j + 1)) 

< jujD-J - ^M\mumD X ju3D]j + ^M\mumD > 

11 1, 1 11 1, 1 

<2^'2^-2l^-2^^><2^'2^^-2l^ + 2^^> 
1/2 ju ~ 1/2 JU \ / 1/2 JU - 1/2 JU 
j 

i2g{2ju + 1) ^ (2jf sin(t52j) 

J = |jC/-iDl 

< juiD\j - lM\mumD >< ju3D;jM\mumD > 

< Ij - \: \m - IM X \j - \ \m - Ib-M > 

1/2 JU - 1/2 1 r 1/2 jt, - 1/2 JU 

3D j - 1 J - 5 J 1 3D j j-\ 

In the last line of Eq. (jBlOjl . we have just done a shift of the summed variable j i.e. j j 



JD 3~\ 3 \ \ 3d 3 + \ 3 } (BIO) 



a. rhs'-^'^') 



^„32 _ 23 n 

rhs^Ng[P / ^23 1^32) (Bll) 
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Using the two-body density matrix Eq. p5|) . 



r/is(5'5) =N J2 ~V^A{J, J'){2jD + 1) < Ij; Om| J'Af > 
J, J' J 

h iu-h ju ) ( h iu - \ iu 

J ^ J I 1 J J' 

+ V3A{J, J'){2jD + 1) < Ij; Oto| JM > 

1 ■ I ■ \ ( 1 ■ 1 ■ 



(B12, 



where 



2 Ji^ 2 ( \ 2 2 

i ^' J I 1 J ^ , 

{1 • 1 ■ 

2 3U - 2 JU 

1 JD-I 3D 

i 3 

A ju^h 3u\ { h ju-h iu 

{A{J,j) { I j^-l \ - AU, J') I I JD-I 3D 

\ 3 J \ 1 J J 



^(J, J') = gexp{-it5(J(J + 1) - J'(./' + l))}(2jy + l)(2j + 1) 
< ju3D;JM\mumD >< 3u3d\ J' M\mumD > 



(B13) 



As J and J' range over the same values, by careful rearrangement of Eq. (|B12p we can reduce rhs^^'^'' to just a 
summation over J, j. 

{h 3u - h 3u ] ( h 3u - \ ju 
I 3D-1 3D } { I 3D - \ 3D 
3 J J I 1 3 J' 

{A{J,3)-A{3,J)) 

= ~NV3gi23D + l){2ju + l)2i ^ sin(tg( J( J + 1) - i(j + l))(2j + 1) ^g^^^ 

JJ 

< 3u3d; JM\mumD >< 3u3d; J'M\mumD >< Ij; Om\J'M > 
k 3u-^ 3u \ ( ^ 3U ~ \ 3U 

2 3D - 2 3D >< 2 3D - 2 3D 

3 3 j [l 3 J 

The 9j-symbol restricts J € {|j — 1|, j + 1}. With these substitutions for J and the help of Eq. (jD8|) and Eq. (|D9|) we 
get, 

r/is(5.5) = i4:V3Ng{2ju + l){2jD + 1) 

X! (2.?' + 1) < jujD;jM\mumD > 

j=\ju~jD\ 

vmTT) { 3D 3D - \ 3] ^^1^) 

{sm{tg2{23 + 1)) < 3u3d\3 + lM\mumD >< Ij; OAf|j + IM > 

{{3 + 1){A'}{B'}{C'} + {3 + i){V'}{£'}{T'}) 
+ sin(-2jtg)) < 3u3d;3 - lM\mumD >< Ij; OM\j ~ IM > 

i3{g}{n}{i} + {3-i){j}{ic}{m 
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where 



^ ju - h ju \ rH'i _ / h JD-h in] trn ^ I 1 J J + 1 



/ /"i . ^ 1 1^ {C} 



J + 1 J + § i ^ \3U-^ 3 + i 3 \ '"'13 + 11 \ 



^ ju - h 3U \ ruT _ / ^ - h 3D \ r-n _ J ^ 3 3 + '^ 



11 1 



2 J I ■J'^ 2 -J ' 2 J J I J ' 2 2 2 



(B16) 



3 1 1 



ic J + 1 j + I J i jc/ - 2 i + i J J I J + 2 2 2 

It can be shown that Eq. (jB15[) is equivalent to, 
rhs^hi) ^ i2V6Ng{2ju + l)(2jD + 1) 

ju+jo 

J2 {^)3u+jo+33+i sm{tg2j) < 3u3D;jM\mumD >< 3u3d]3 - lM\mumD > 

3 = \ju-jD\ 

VW^l { ^""jD ^ 3D- i J ^ 1 } < " ^' ^^^^^^^ > 
(jM}{6}{C}+(j + l){I?}{£}m) 

(j{g}{7i}{x} + (j-i)m{/c}{/:}) 

where, 



3u - 2 ju 1 TRi J 2 - 2 I sr\-] 1 3-^3 



^ J(7 - ^ J't' \ rci _ / ^ ic - ^ 1 r-r-1 _ / 1 J'-I J 



(BI8) 



The important thing about Eq. (jB17p . is that it has exactly the same summation limits as Ihs^^'D Eq. (jBlOp . So if 
we define, 

ju+jo 
j=\3u-jD\ 

ju+jo 

r/is(5.i) = ^ Y (B20) 

j=\ju-jD\ 

then with the help of Eq. jDlJ-Eq. (jD4|) and after some simple but tedious algebra, one arrives at, 

(_)-2(2,+,.+,.)((_)2, _ (_)2M)(^. 



X -Y = 



2j3u{23u + l) 

I (B21) 

{3 + 3D - 3U){j - jo + jU){l - j + JD + 3UKI + 3 + 3D + 3U){3 - M){j - M) 



Y (2j-l)(2j + l) 

Now if we impose the physical constraints. 



3 =3u-Jd+c (B22) 
M = Ju - 3D (B23) 



23 



where c is a positive integer, then, 

= (_)20£/-Jd)^_)2c _ ^-^ (B24) 
= 

Hence, 

X -Y = (B25) 

Thus wc have shown the final piece of the proof, that Ihs^^'^^ = r/is^S'l). Therefore we have proven the consistency 
of the density matrix with the kinetic equation. 
Q.E.D. 



APPENDIX C: ONE AND TWO BODY DENSITY MATRICES WITH INTERACTION AND VACUUM 

OSCILLATION 



Pi it) = E (-)■'"■'' exp{-ii[B(M^ - M^') + g{J{J + l) - J'( J' + l))]} 

J,J',M^,M^' 

B B B' 

<B,M(-/3)<,MB'(-/3)C„.f (-/3)Cf' ,,r.,'(/5')(2ic/ + l)(2i + 1) 

< jujo; JM\mumD >< jujo; J'M\mumD > (CI) 

< iii; JM^Imf X jij; J'M^|mf' > 

/ Jl ic/ " h 3N \ / Jl ju - jl JN \ 
\jD J J j\jD J' j j 

liimi >< jim[\ 

Pi2(i)= E exp{-z<[i3(M^-M^')+5(J(J + l)-J'(J' + l))]} 

j,j',mb,mB' 

ii2,Ji2',j,m^ 

™?2>™12 i'"12,"ll2' 

(2it, + + l)(2i + l)[(2ii2 + l){2j[^ + 

< jujo; JM\mumD >< jujo; J'M\mumD > 

< ii2j;mf2m^| JM^ >< j[^j-mf^m'^\J'M'^' > 

jl ku ju 1 [ jl ku ju 1 
j2 ko jo > I j2 ko jo \ 

ii2 j J ] \ ii2 3 J' J 

|jl2mi2 ><ii2"^'l2l 

/3 is the angle of rotation between the flavour axis and B. Spin projection superscripted with B is the projection 
in the rotated basis. B is the magnitude of B. d'j^B m ^ rotation matrix. All other parameters are as defined in 
§sec:summaryDef. 



APPENDIX D: ANGULAR MOMENTUM FORMULAE 

3j-symbol formulae: 
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fa b a + b\ ^r.a-b--, I (2a)!(2b)!(fl + 6 + 7)!(a + 

\a S 7 J ^ > y {2a + 2b + iy.{a + ay.{a ^ ay.{b + (3y.{b - py. 



(D2) 



a 6 7 ''^^^^ "^^Y (2a + 26)!(a + a)!(a-a)!(6 + /3)!(6-/3)! ^^^^ 



a b a + \ _ , ^a-h-7o/.„ / (2a - 1)!(26 - l)!(a + 6 + 7 - l)!(a + 6 - 7 - 1)! 



a 6 a + & - 2 \ , xa-b-^ / + ^ - 7 - 2)!(a + b + 7 - 2)!(2a - 2)!(2b - 2)! 

a <5 7 y ^ V 2(a-a)!(a + a)!(5 + /3)!(6-/3)!(2a + 26-l)! ' 



6j-symbol formulae: 



(a6)e,d;c>=^|a,(M)/;c>|(2e + l)(2/ + l)]^H''+''+^+'^|^ J ^ 



a b e 
d c f 



9j-symbol formulae: 



a b c 

d e c = 



g g l\ V4c(c+l)(2c+l)5(ff+l)(2ff+l) 
These angular momentum formulae are from 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7 

[s: 

[9; 
[lo; 
[11 

[12 



(D5) 



{-Y+''+''+''W{abcd;ef) (D6) 



a b c 

< {ab)c, {de)f; i\iad)g, {be)h- 1 >= [(2c + 1)(2/ + 1)(2.9 + l){2h + Ip { d e f } (D7) 

g h i 



d e f K^^ ^s'^^ > W{abcd;ef) (D8) 

„ V(2c+l)(2ff+l) 



+ - + ^) - + ^) + '^(^ + ^) (-)(-+^--) W^(abcrf; e/) (D9) 
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